Infrared magneto-spectroscopy of graphite in tilted fields 
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The electronic structure of Bernal-stacked graphite subject to tilted magnetic fields has been inves- 
tigated using infrared magneto-transmission experiments. With the increasing in-plane component 
of the magnetic field Bu , we observe significant broadening and partially also splitting of interband 
inter-Landau level transitions, which originate at the H point of the graphite Brillouin zone, where 
the charge carriers behave as massless Dirac fermions. The observed behavior is attributed to the 
lifting of the twofold degeneracy of Landau levels at the H point - a degeneracy which in graphite 
complements the standard spin and valley degeneracies typical of graphene. 

PACS numbers: 71.20.-b, 71.70.Di 



I. INTRODUCTION 

It was the fabrication of single-layer graphene 1 -^ and 
subsequent discovery of massless Dirac fermions 3 -! 4 - which 
triggered the present increased interest in the electrical 
and optical properties of graphite 5 - - — - supposedly a well- 
known material for the condensed matter physics. 

Even though graphene is a purely two-dimensional 
(2D) system and graphite is characterized by a (highly 
anisotropic but still) clearly 3D band structure, these 
materials, as demonstrated experimentally^ ' 10 ' 18 share 
surprisingly similar optical response when the magnetic 
field B is applied perpendicularly to layers. A sim- 
ple model, invoking inter-Landau level excitations be- 
tween highly-degenerate Landau levels (LLs) of mass- 
less Dirac fermions, implies the magneto-optical response 
that is linear in y/~B, see, e.g., Refs. [l9l - |26l and is ca- 
pable to account for a significant part of the magneto- 
optical data acquired on graphite. Importantly, these 
data come not only from recent magneto-transmission 
studies of thin specimens^ ' 10 ' 11 but also from original 
measurements carried out in late seventies^ in which 
\/S-scalcd spectral features have been observed using the 
magneto-reflection technique. This pioneering work is a 
good candidate for the first direct experimental observa- 
tion of massless Dirac fermions, which in bulk graphite 
coexist with massive particles and which provide more 
conventional, i.e., linear in B response ! 28 ' 29 

The electronic band structure of graphite in the 
magnetic field is mostly described using the standard 
model proposed by Slonczcwski, Weiss and McClurc 
(SWM) ) 30 ' 31 even though presumably more precise, but 
at the same time, also more time-consum ing approaches 
appeared recently, see, e.g., Refs. l32l andl33l The SWM 
model has been derived in late fifties using mostly sym- 
metry arguments; it describes the electronic structure 
near the H-K-H edge of the Brillouin zone with energies 
not too distant from the Fermi level. The six of seven pa- 
rameters in the SWM model, 70, • ■ • , 7sr^ are usually in- 



terpreted as tight-binding hopping integrals between the 
nearest-neighbor and partially also next-nearest-neighbor 
atoms. An additional parameter A, related to the non- 
equivalence of carbon atoms in A and B positions, is re- 
ferred to as a pseudogap. All parameters must be consid- 
ered rather as adjustable parameters than true hopping 
integrals and arc usually obtained by fitting either exper- 
imental data or the results of ab initio calculations^ The 
importance of individual parameters significantly varies, 
depending on the type of experimental data for interpre- 
tation of which the SWM model is used. 

For instance, the periods of Shubnikov-de Haa o 6 ' 36 " — 
(SdH) and de Haas-van Alphen^ - (dHvA) oscillations 
depend on the extremal cross sections of the complex 
Fermi surface and all SWM tight-binding parameters 
must be properly taken into account. Similarly, cyclotron 
resonance experiments , 28 ' 41 ' 42 which are also sensitive to 
the immediate vicinity of the Fermi level and which pro- 
vide fairly rich response, can be hardly understood with- 
out the full SWM model. 

On the other hand, interband transitions between elec- 
tronic states far away from the Fermi surface can be 
successfully described using a simplified approach ; 10 i 11 
which models the magneto-optical response of bulk 
graphite as a sum of responses of an effective graphene 
bilayer and monolayer^ 4 - Notably, the physical prop- 
erties of a 3D system are thus described using re- 
sponses of two purely 2D materials, and interestingly, 
not more than two coupling constants, intralayer 70 and 
interlayer 71, are needed in the very first approach.4 5 - 
Within such a minimal model, the H point provides 
response similar to a single graphene sheet, but richer 
due to an additional twofold degeneracy, and the K 
point behaves as bilayer graphene, however, with the 
interlayer coupling enhanced twice as compared to the 
true bilayer. Limits of this model have been found, 
e.g., by revealing the electron- hole asymmetry at the K 
point of bulk graphite in recent magneto-transmissionji 1 - 
magneto-reflectio n 14 ' 17 and magneto-Raman studies. 15 ' 16 
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FIG. 1. Transmission of a thin graphite specimen measured at magnetic field tilted by tp — 0°, 15° and 30° with respect to 
the c-axis of graphite. Values nearby individual curves always denote component of the magnetic field perpendicular to the 
graphite layer B± . The missing parts of spectra correspond to regions in which the supporting tape is completely opaque. The 
energy axis for each curve is rescaled by a factor of \/B± to facilitate identification of spectral features originating at the H 
point of graphite. Vertical lines correspond to positions of van Hove singularities in the joint density of states as calculated 
using the minimal nearest-neighbor tight-binding model.— Since the lifting of degeneracy exactly at the H point is governed 
by the coupling constants T n , n +i oc \ / B± tan(y>), the positions of these van Hove singularities remain constant in these figures. 
For clarity, successive spectra in (a),(b), and (c) are shifted vertically by 0.1. 



The full SWM model has to be used in such a case to get 
quantitative agreement between the experimental data 
and theory. 

In this paper, we set other limitations of the effective 
monolayer and bilayer model for the magneto-optical re- 
sponse of graphite. Namely, we test its applicability in 
experiments performed in the tiltcd-ficld configuration, 
B = (0, B|| , B±), which is a basic tool to distinguish be- 
tween 2D and 3D character of condensed matter systems. 
We focus on the graphene-like signal from the H point 
and show that the magneto-optical response of graphite 
in a tilted magnetic field follows the total magnetic field 
and not only its perpendicular component, as should be 
in the case of an ideal 2D system. The infrared magneto- 
transmission technique is thus, perhaps surprisingly, sig- 
nificantly more sensitive to the in-plane component of the 
magnetic field as compared to other techniques such as 
SdH or dHvA oscillations, which reveal the 3D charac- 
ter of graphite only for rather high tilting anglesi 38 ' 40 To 
interpret our data, we use recently developed theory of 
the graphite band structure subject to a tilted magnetic 
field, which predicts lifting of the twofold degeneracy at 
the H point M- This degeneracy, taking origin in the 3D 
character of graphite (four atoms in a unit cell instead of 
two for graphene), is an additional one to the valley and 



spin degeneracies in graphene. 



II. EXPERIMENT 

Thin graphite specimens for our magneto-transmission 
study have been prepared by exfoliation. A thin layer 
of bulk graphite, with an average thickness around 
wlOO nm, was located on the scotch tape used for ex- 
foliation, which has several relatively wide spectral win- 
dows with a sufficiently high optical transmission. A 
high-quality natural graphite crystal has been chosen 
for exfoliation, since it provides equivalent but better 
pronounced magneto-optical response as compared to, 
e.g., highly-oriented pyrolytic graphite^ The magneto- 
transmission spectra were measured for the magnetic 
field inclined with respect to the c-axis of graphite by 
selected angles of (p, i.e., in the perpendicular (ip = 0°) 
and several tilted-field configurations. 

To measure the transmission spectra in the spectral 
range 100-800 meV, the non-polarized radiation of a glo- 
bar was analyzed by a Fourier transform spectrometer 
and guided to the sample by light-pipe optics. The sam- 
ple was placed in a cryostat at temperature of 2 K lo- 
cated inside superconducting and resistive coils, which 
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FIG. 2. Relative magneto-transmission spectra for tilting an- 
gles tp = 0°, 15°, 20°, 25°, 35° and 50° taken with the same 
perpendicular component of the field B± = 7 T. The broad- 
ening of the lines with increasing angle, i.e., with the increas- 
ing in-plane component of the field Bit is well visible, e.g., 
on the C line, the width of which increases roughly linearly 
with B«. A sample different from Fig. [T] has been used (with 
higher density of crystallites and their average thickness). For 
clarity, successive spectra are shifted vertically by 0.1. 



reach magnetic fields up to 13 T and 28 T, respectively. 
The transmitted signal was detected by the composite Si 
bolometer. All spectra presented in this study have been 
normalized by the zero-field transmission. 



III. RESULTS 

The magneto-transmission spectra measured at three 
different angles between the c-axis of graphite and mag- 
netic field, tp = 0°, 15° and 30°, are presented in Figs.QJ,- 
c, respectively. To facilitate identification of individual 
absorption lines, the transmission curves are depicted as 
function of the photon energy normalized by the factor 
of y/B±, which is typical of LLs in a system of ideal 2D 
massless Dirac fcrmions. Plotted this way, we can iden- 
tify graphene-like signal originating at the H point and 
easily follow its deviation from the y/B± dependence in- 
duced by the in-plane component of the magnetic field 
J3||. An additional set of data is presented in Fig. [5J 
where the magneto-transmission spectra for several tilt- 
ing angles tp are presented, all measured at a fixed per- 
pendicular component of the field = 7 T. 

The observed absorption lines have been marked con- 



FIG. 3. Landau subbands in the vicinity of the H points, 
k z d = 7r/2, at B± — 5 T. The solid and dashed lines describe 
the k z dependence of energies E„ and E„ in even and odd 
graphene sheets at (p = 0°, 15° and 30°. 



sistcntly with the notation introduced earlier.— ^ The 
transitions denoted by Roman letters have their direct 
counterpart in the response of graphene^ while the 
"Greek" lines are characteristic of graphite.— They are, in 
principle, dipole-forbidden in a pure 2D system of Dirac 
fermions, nevertheless, they can be consistently explained 
with the standard dipole selection rule (n — > n± 1) when 
the twofold degeneracy of LLs at the H points of bulk 
graphite is properly considered. Let us note that the 
widths of absorption lines reflect not only the naturally 
present disorder in the graphite crystal, but they are also 
partially given by the profile of the individual Landau 
subbands in the vicinity of the H points 

The transmission spectra measured at ip = 0°, pre- 
sented in Fig. Q}t, are fully analogous to previous 
experiments^ ' 10 ' 11 and the color vertical lines mark trans- 
mission minima proportional to y/B±. All such transi- 
tions originate at the H point. The transmission minima 
are more pronounced at lower energies and their width 
increases with the increasing magnetic field. Interest- 
ingly, the number of observed transitions remains nearly 
constant with B± . This behavior reflects the specific en- 
ergy dependence of the relaxation rate (i.e., broadening 
of lines), t~ 1 (E) oc \E\, which maps the (linear in energy) 
density of states around the H point. An analogous effect 
has been recently observed also in graphene specimens.— 

At a non-zero tilting angle tp, the observed magneto- 
transmission spectra significantly deviate from expecta- 
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tions for a purely 2D system, which is in the case of 
orbital effects only sensitive to the perpendicular compo- 
nent of the field. At ip as low as 15°, sec Fig.QJi, the tran- 
sitions denoted by Roman letters change the shape and 
broaden, while the "Greek" lines become significantly 
weaker. For ip > 30°, the "Greek" lines completely dis- 
appear from spectra and also the Roman lines are much 
less pronounced as a result of a significant broadening. 
Alternatively, we can follow these effects in Fig. [2j where 
the magneto-transmission spectra are plotted at several 
angles with the perpendicular component of the field kept 
constant, B± = 7 T. The effects induced by the in-plane 
component of the field Bu are well illustrated, e.g., on 
the C line. With increasing tilting angle, this line does 
not only significantly broaden, but also gains a complex 
structure - a strong asymmetry is developed and the line 
becomes nearly split into two components for higher an- 
gles (tp > 25°). 

IV. DISCUSSION 

To interpret the broadening of absorption lines with 
Bii, we will consider the electronic band structure at the 
H point of graphite in detail. In particular, we will fo- 
cus on the Bt i -induced lifting of the twofold degeneracy, 
which in graphite complements the spin and valley de- 
generacies in graphene, and follow the theory recently 
developed by Goncharuk and Smrcka4^ The additional 
twofold degeneracy may be interpreted as a direct con- 
sequence of the effectively vanishing interlayer interac- 
tion for a charge carrier with the momentum k z d = n/2, 
where d is the interlayer distance. The reason is that 
the neighboring graphene sheets are rotated by 30°. If 
the field dependence of the energy in even layers is 
E e n oc ±\/B±n then the energy in odd layers reads 
E° oc ±^Bj_(n + 1), where n = 0,1, - ■ ■ is the index 
of the LLs. Two states \n + l) e and \n)° belonging to de- 
generated eigenenergies E^ +1 = E° are orthogonal and, 
therefore, the corresponding interlayer hopping integral is 
equal to zero for B\\ = 0. In tilted fields the in-plane field 
component B\ \ shifts the mutual position of centers of or- 
bits |n + l) e and \n)° in real space by yo = dBu/B±. The 
orbits are no longer exactly orthogonal and the interlayer 
interaction does not completely vanish. In the lowest or- 
der of the perturbation theory we get instead of E° and 
E° n+1 four energies E^' = ±T n , n+1 ± [(£°) 2 +T„ 2 „ +1 ]V 2 

where T n , n+1 = 71 d tan(</>) y/B ± \e\{n + 1)/{2K}M- Let 
us note that another theory presented in Ref. |47| is de- 
voted to the case of magnetic field applied strictly parallel 
to the sheets of the graphene bilayer and multilayers in- 
cluding graphite. It is suggested that the obtained energy 
spectrum can be verified experimentally using electron 
tunneling or optical spectroscopy. 

Obviously, the first order perturbation theory em- 
ployed in Ref. HH, which involves only two states with the 
same energy exactly at k z d = 7r/2, is not the best approx- 
imation and overestimates the splitting. It is acceptable 



only for rather small tilting angles and LLs with limited 
n. For this reason, we have calculated the eigenenergies 
numerically employing a larger basis, which allows us to 
calculate also the k z dependence in the vicinity of the H 
points. The results are shown in Fig. [3] Dashed lines de- 
scribe the k z dependence of energies E^ and E° in even 
and odd graphene sheets at ip — 0°. Solid lines describe 
the energies E^ and E° , both at ip = 15° and 30°, re- 
solved by different colors. The corresponding cigenstates 
mix the wave functions from both even and odd layers, 
nevertheless, the same notation is kept to emphasize to 
what energies E^ and E° are reduced for ip — > 0. 

Our numerical calculations reveal the importance of 
the states in the vicinity of the H point. Two side ex- 
trema in the k z dependence of energy subbands appear 
and the curvature of the k z curves is reversed at the H 
point. Each new local extremum developed along the k z 
dependence of the energy subbands adds a new van Hove 
singularity to the (joint) density of states. The energy 
gap opened between subbands appears slightly away of 
the H point and it decreases with the increasing LL index 
n (n > 2). Unlike the splitting directly at the H point, 
this energy gap of is no longer proportional to \/B±. 

At ip = 0, each absorption line consists of up to four 
degenerate transitions. To be more specific, we have two 
transitions for each "Greek" line and also for the B line, 
the rest of lines denoted by Roman letters include four 
degenerate transitions, with corresponding four (degener- 
ate) van Hove singularities in the joint density of states, 
for details see Refs. H and At <p = 0, positions of 
these singularities are represented in Fig. QJi by the ver- 
tical color lines. The situation becomes more complex 
at tp > 0, when the degeneracy of these four Van Hove 
singularities is lifted. In addition, other singularities, pre- 
sumably weaker, develop in the vicinity of H point. To 
illustrate the strength of effects induced by the in-plane 
magnetic field, we plot, in Figs.[Tb,c, positions four main 
singularities in the limits of validity of the first order ap- 
proximation - those which originate directly at the H 
point. Two of them remain almost degenerate. On the 
other hand, for the sake of clarity, we do not mark po- 
sitions of additional singularities developed due to anti- 
crossing of Landau subbands further from the H point. 
To justify this, we note that such singularities in princi- 
ple exist even at B\\ =0 due to the trigonal warping term 
which is neglected in the simple effective monolayer 
and bilayer model. Nevertheless, they have not been ob- 
served in the experiment in perpendicular fields.— 

To sum up, the in-plane magnetic field profoundly 
modifies the profile of Landau subbands at the H point 
of bulk graphite. A scries of minigaps is created al- 
ready within the first order of the perturbation theory, 
which is directly reflected by newly developed van Hove 
singularities in the joint density of states. Experimen- 
tally, this leads to the splitting of the observed dipole- 
allowed transitions with the increasing angle tp, or at 
least, if the disorder effects are realistically involved, to 
a significant broadening of these transitions. This split- 
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ting/broadening increases at the fixed B± roughly lin- 
early with By , or equivalently, approximatively scales as 
\/B± when the tilting angle if is kept constant. Such be- 
havior is consistent with the experimental data plotted 
in Fig. [5] and in Figs.QJi-c, respectively. 

V. CONCLUSIONS 

The electronic band structure at the H point of bulk 
graphite has been studied using the infrared transmis- 
sion technique in magnetic fields tilted with respect to 
the c-axis of this material. While for the magnetic field 
applied along this axis, the magneto-optical responses of 
graphite (due to the H point) and graphene (due to the 
K point) closely resemble each other, pronounced devi- 
ations clearly appear with the increasing tilting angle. 



The 3D nature of the electronic band structure of bulk 
graphite is thus revealed at significantly lower tilting an- 
gles as compared to SdH and dHvA measurements! 38 ' 40 
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